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One reports computational study revealing a set of general requirements, fulfilling of which would 
allow employing changes in ambient conditions to regulate accomplishing the collective outcome of 
emerging active network patterns in an ensemble of a moderate number of finite discrete systems. 
The patterns within all these component systems emerge out of random expansion process governed 
by certain local rule. The systems modeled are of the same type but different in details, finite discrete 
spatial domains of the expansion within the systems are equivalent regular hexagonal arrays. The 
way in which elements of a component system function in the local information transmission allows 
dividing them into two classes. One class is represented by zero- dimensional entities coupled into 
pairs identified at the array sites being nearest neighbors. The pairs preserve their orientation 
in the space while experiencing conditional hopping to positions close by and transferring certain 
information portions. Messenger particles hopping to signal the pairs for the conditional jumping 
constitute the other class. Contribution from the hopping pairs results in finite size effects being 
specific feature of accomplishing the mean expected network pattern representing the collective 
outcome. It is shown how manifestations of the finite size effects allow using changes in parameters 
of the model ambient conditions of the ensemble evolution to regulate accomplishing the collective 
outcome representation. 
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I. INTRODUCTION 

This study contributes to a search for methods of reg- 
ulating the collective outcome of assembly processes tak- 
ing place in an ensemble of a moderate number M of 
finite complex systems, of the same type but different in 
details, evolving simultaneously and independently, un- 
der influence of stochastic factors within equivalent finite 
spatial domains and in the same ambient conditions. The 
finite discrete spatial domains are equivalent regular ar- 
rays in all the component systems. Information transmis- 
sion processes within these domains result in manifesta- 
tions represented by various patterns of expanding spatial 
regions specified with the same feature in all the systems 
at each stage of the process. One assumes that these re- 
gions specified within all the component systems affect 
simultaneously certain receiver system. Moreover, differ- 
ences between contributions, which distinct component 
systems have to the effect of their influence on the re- 
ceiver system, are attributed only to differences between 
spatial arrangements of the regions specified within these 
component systems. Accordingly, spatial arrangements 
of these regions within all the systems correspond to the 
receiver system transformation that manifests itself as a 
collective effect (CE) at each stage of the process. Thus 
the receiver system recognizes collective outcome of evo- 
lution processes within all the component systems. One 
also assumes that all the component systems contribute 
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to the collective effect in the equivalent conditions. For 
the reasons just noted, the CE can be represented by 
certain pattern of distinguished regions within a spatial 
domain being equivalent to the domain of a component 
system. Changing certain values parameterizing the am- 
bient conditions of evolution of all the component sys- 
tems is the way in which regulation of the collective effect 
is under consideration here. 

In course of this study we have found that a prereq- 
uisite for following this way is presence of finite-size ef- 
fects (we use abbreviation, F-SE) i.e. specific features of 
the CE accomplishment which result from realizing the 
local information transfer in the component systems of 
the evolving ensemble by one dimensional, finite-size el- 
ements able to change their positions in finite discrete 
space while transmitting certain signals between their 
ends and preserving their orientation in the space. The 
signal is considered as finite portion of information trans- 
mitted discretely between an element ends. Accordingly, 
each such element is represented by the pair of its ends 
only, one end that only sends a signal and the other re- 
ceiving it. These ends coincide with sites of the regular 
array which are nearest neighbors. Hopping of the pair 
to a separate position close by requires obeying certain 
stochastic condition. Whole the pattern within each com- 
ponent system is being accomplished at each step of the 
evolution and can be composed only of sites being desti- 
nations for the receiving ends of the pairs that have expe- 
rienced hopping at this step. This results in high diver- 
sity of the patterns which manifests itself in appearance 
of rugged islands in course of a component system evo- 
lution initiated from a single cluster; this specific diver- 
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sification is a prerequisite for appearance of the F-SE. It 
has been also found that manifestations of the F-SE allow 
regulating the CE accomplishment effectively, in the way 
just mentioned, if elements constituting other class signal 
the pairs to perform the conditional hopping. This class 
encompasses elements represented by zero-dimensional 
entities transferring the information portions by hopping 
between closest neighbor sites. The way of modeling the 
local information transfer just mentioned is explained in 
Section [n]and algorithms of modeling the conditional in- 
formation transfer are presented in Section IIIII This has 
been assumed as simplest model allowing the manifesta- 
tions of effects being the subject investigated in course of 
this research (this is elucidated in Sections IIVI and IVIII . 

Remarkable manifestations of the finite-size effects, 
finiteness and discreteness have appeared indispens- 
able to allow regulating the CE accomplishment. The 
number, M of the component systems affects usefulness 
of these manifestations for the regulation. The term mod- 
erate M used hereafter means: M is so large that the col- 
lective effect is accomplished reliably with certain specific 
features and M is so small that all diverse specific man- 
ifestations attributed to the discreteness, finiteness and 
transfers resulting in the finite-size effects are remarkable 
and can be easily identified in the CE. The way of esti- 
mating the lower bound Mq of the M being moderate 
and role of other values parameterizing the model am- 
bient conditions are presented in Section IIVI Criterion 
of detecting the upper bound, Mi to the value M con- 
sidered as being moderate and its role are explained in 
Sections EUlEini)- 

One of the motives for doing this research is the ne- 
cessity in predicting the opportunity for regulating the 
collective effect emerging out of evolution of an ensem- 
ble of systems of nanoelements that are assembling quasi 
two-dimensional functional networks over prepatterned 
surfaces (for the purpose of modeling this process, the 
nanoelements can be represented by pairs of their ends 
only in the way alluded to previously). This task may 
be considered as being relevant to low cost process of 
self- assembling nanostructured functional materials (this 
is the so called bottom-up method) that are expected to 
be applied as sensors, electric field emitters stimulating 
living tissues and as other active covers. Then, an elec- 
tron would be usually the signal transmitted by the na- 
noelements and this explains the convention of identi- 
fying them as nanowires. The short mobile nanowires 
would differ from one another in length and conductiv- 
ity, and the conduction would influence binding other 
elements at their ends. A number of them would be teth- 
ered, either temporarily or permanently, by elements of 
the surface texture that would constrain freedom of their 
displacements and the differences between the nano- wires 
would contribute to the system randomness. One may 
suggest a number of candidates for the material com- 
ponents of such systems. For example, properties of 
single-walled carbon nanotube conductors (sec also 
http://cnanotech.com) allow considering their short 



intervals as candidates for the short nanowires. The 
short nanotube intervals can emit electric field from their 
ends, form stable colloidal suspension in water Hand can 
be adsorbed spontaneously to a gold surface 0]. Pos- 
sibility of fabricating multivalent receptors on the gold 
nanoparticle scaffolds Q seems to open the opportunity 
for employing them to recognize elements self-assembling 
into the network. Moreover, gold nanocrystals can self- 
assemble into ordered hexagonal texture of the surface 

and electron transfer appears to be engaged into their 
interactions with fullerene particles Q . This reveals fea- 
sibility of a material system evolution being relevant to 
the model process reported here. Let us note, however, 
that general character of the model assumptions made 
in this research may make the presented approach useful 
for studying other physical processes also. 

A computational study of the opportunity for regu- 
lating the collective effect emerging out of the ensemble 
of evolving complex systems can be carried on by sim- 
ulating the model process proposed here. We employ 
the known idea Q, 0, U E3 of considering signals being 
finite portions of information about certain generalized 
model feature transported in discrete way between sites 
of a finite discrete space (the space is represented here 
by finite regular hexagonal lattice). Appearance of that 
model feature at a site is indicated by covering this site. 
Thus a pattern of covered sites within spatial domain of 
a system represents effect of the system evolution. 

It is specific to this model that the scheme of cover- 
ing the sites is assumed so as to represent effectively the 
underlying transport mechanisms resulting, eventually, 
in the finite-size effects. For this purpose we adapt our 
method of discrete displacements, within frame of which 
hopping pairs of zero-dimensional entities have been em- 
ployed to simulate effects of information transmission rel- 
evant to turbulent transport and, then, agreement with 
available experimental results has been demonstrated [j| . 
Here, the collective effect is represented by the pattern 
being mean expected form of the patterns emerging un- 
der equivalent conditions, simultaneously and indepen- 
dently (in parallel) in all systems constituting the en- 
semble. 

In this work, one investigates random expansion pro- 
cess (REP) as Markov process of covering sites of finite 
regular hexagonal array. States of the REP at stages 
T are finite random sets (FRS), F(T) = {Ai{T),i = 
1,2, ...,M}. The realizations Ai(T) develop simultane- 
ously and independently in the same model ambient con- 
ditions. Accomplishing of the model collective effect is 
represented by sequence of patterns J F(T), each of which 
is mean expected form of all the Aj(T) and computed 
as mean measure set (MMS) after known algorithm [Tl| 
briefly reported in the Appendix D 

Characteristics of the MMS evolution pattern are pre- 
sented in Section Specific, predictable changes in the 
MMS evolution patterns, which result from varying cer- 
tain parameters of the model ambient conditions, have 
been recognized as features attributed to the finite-size 
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effects modeling. They reveal a possible way of searching 
for a method of regulating the collective effect emerging 
out of a material system ensemble for which the process 
underlying evolution of the ensemble may be adequately 
represented by the REP. These results are presented and 
discussed in Section IVll 



II. METHOD OF MODELING THE LOCAL 
INFORMATION TRANSFER 

The way of modeling the mechanism of information 
transfer in each realization Ai evolving from a stage T 
to T + 1 can be elucidated by considering certain ideal- 
ized physical system. Let us assume that the system 
consists of short nanowires, transferring signals while 
hopping to positions close by, and of certain complex 
molecules occupying all sites of the finite hexagonal ar- 
ray representing texture of a prepatterned surface. More- 
over, transmission of the signal changes binding proper- 
ties at the nanowire ends and this can affect behavior 
of the complex molecules occupying sites coinciding with 
these ends. With the purpose to model effects produced 
by the nanowires, we generalize our approach used 
previously for modeling the contribution which certain 
finite-size elements representing organized fluid streams 
and structures have to the collective effect of transport 
processes in the turbulent shear flow. In accordance with 
this approach, employed to situation being considered 
here, the finite-size elements transferring the signals are 
modeled as pairs of neighbor sites (x, y) that transmit 
a signal from x to y while hopping effectively to a po- 
sition (x',y') close by. This effective hoppings are real- 
ized so that information portion about orientation of two 
neighbor sites coupled into a pair is transfered from the 
position (x, y) to (x',j/) (see the triplet of patterns in 
the bottom of Fig.^Jt). Generally, each pair of neighbor 
sites can be considered coinciding with both ends of a 
model nanowire. Let us explain in this connection that 
the model being presented should not be considered liter- 
ally: The model is effective in character , i.e. local model 
processes can represent cumulative outcome of all pos- 
sible local actual realizations allowed in certain ambient 
conditions by the space configuration (here, hexagonal 
array) and related to one step of the model evolution. 
Further we explain spatiotemporal aspect of the effective 
modeling and thus reveal way of representing diversity 
of possible events which results from participation of the 
hopping pairs in the information transmission. 

The pattern identifying the realization Ai(T) is com- 
posed of sites covered at the stage T only. The covering 
of a site represents activation of a complex molecule at 
the site in the idealized physical system considered. In 
this system, each molecule activated sends identical mes- 
senger particles to all its nearest neighbors. In the com- 
putational model, each messenger particle is represented 
by a signal C being the same for each site covered from 
which the C is sent. Accordingly, each site occupied by 



such activated molecule is identified as a site with dis- 
tributable feature C . One assumes that each molecule 
being active at a stage T sends unconditionally the mes- 
senger particles C to all its nearest neighbors at the step 
of evolution, from T to T + 1, and thus loses its activity. 
A molecule can become active at the stage T + 1 only in 
result of a process identified hereafter as transferring the 
molecule activity by a nanowire at the step of evolution 
from T to T + 1; this process is elucidated below. For 
first, we explain spatial aspect of the effective modeling 
of the local information transfer and this will help us to 
explain subsequently its temporal aspect. 

One assumes that the messenger particle C coming to 
a molecule (either activated or not) from a nearest neigh- 
bor one can be effectively redirected to any molecule be- 
ing the nearest neighbor of this intermediate one (see 
Fig. nji). In virtue of the effective modeling alluded 
to previously, we represent contribution that these local 
redirection events have to effects of the transport pro- 
cess diversity as local cumulative result of the all possi- 
ble model events of the redirection. Note that differences 
between the model events contributing to the local cumu- 
lative result concern only ways of the component redirec- 
tions and are attributed to the space topology identified 
as regular hexagonal tiling. As no direction in this reg- 
ular array is distinguished, the local cumulative result 
can be assigned to a virtual situation in which an active 
molecule would send for six messenger particles to each 
its nearest neighbor molecule that would redirect them, 
for one, to each of its six nearest neighbor molecules. The 
redirection does not cause any changes of the molecule 
redirecting the particles but the messenger particles redi- 
rected acquire new feature allowing them activating each 
nanowire whose sending end coincides with a site occu- 
pied by a molecule being nearest neighbor of the one that 
has redirected the particle. This new feature is consid- 
ered as enhancement of the signal C due to the redi- 
rection and we denote the signal enhanced as C r . Note 
that arrangement of sites into the hexagonal array al- 
lows activating nanowires at the six initial positions by 
one messenger particle C r achieving the site coinciding 
with their overlapping sending ends. Here, for clarity, we 
continue the explanation for example of possible behav- 
ior of a nanowire at one initial position (sufficient condi- 
tions to be obeyed for activating those nanowires as well 
as nature of the effective modeling the transport process 
that engages them are revealed hereafter). A nanowire 
activated by the signal C r delivered to its sending end 
can act as a carrier hopping to a position close by while 
transmitting certain signal Q emitted from its sending 
end toward the receiving one (see Fig. QJi) . The C r does 
not affect the complex molecule at the site coinciding 
with the sending end of the nanowire. Moreover, state of 
this molecule, which can be intact or activated, has no 
effect on activating the nanowire by the C r . The Q itself 
is not a molecule activation signal, transmission of the Q 
by the nanowire is required to change binding properties 
at the nanowire ends: then the sending end can recognize 
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a molecule that was active at a stage T and this allows 
reception of the nanowire at its receiving end by the in- 
tact molecule occupying the same site as this end at the 
stage T+ 1. This reception makes this molecule active at 
the stage T+ 1. The reception of a single active nanowire 
is sufficient to activate the molecule. The molecule that 
can be activated in the way just presented will be identi- 
fied hereafter as target molecule occupying the target site. 
In accordance with assumptions alluded to previously, an 
active nanowire can experience hopping only to a posi- 
tion identified by a pair of neighbor sites so orientated 
in space as the pair of neighbor sites representing the 
nanowire at its initial position. Moreover, one assumes 
that directions of the nanowire displacements agree with 
direction of transmitting the Q, from the nanowire send- 
ing end to the receiving one. Accordingly, there are five 
such positions close by which are identified by pairs of 
neighbor sites that do not coincide with any of the two 
sites identifying the initial position (see Fig.^Jt). Let us 
recall that, within frame of the information transfer rep- 
resentation being considered, the state of transmitting 
the signal Q and hopping of the active nanowire take 
place at one step of evolution which is elementary unit of 
time here. Therefore, change of binding properties at the 
active nanowire ends due to appearance of the nanowire 
in the state of transmitting the signal Q concerns also 
the nanowire at its initial position. The nanowire action 
following the change in binding properties of its ends can 
be considered as direct transfer of the molecule activ- 
ity. Eventually, we consider an active nanowire that can 
transfer directly the molecule activity at any of the six 
positions identified as destination positions. Thus it can 
contribute to the information transmission stream that 
is initiated by emitting the messenger particles, C and 
occurs at each step of the evolution. 

Let us explain in this connection the nature of the 
effective modeling of the local transfer realized by the 
nanowires. Our approach is like the one used in the case 
of the messenger particle redistribution elucidated here 
previously. Contribution that the local hoppings of the 
active nanowires have to effects of spatial diversity in the 
transfer processes for one step of the evolution, from the 
T to T + 1, is modeled as result of a virtual event in 
course of which six nanowires are activated at one ini- 
tial position, experience hopping and are able to transfer 
directly the molecule activity at the six destination po- 
sitions respectively (see Fig. [TJ,) . Let us note that the 
pattern just sketched for example of a nanowire activated 
at one initial position remains the same for each of six 
possible initial positions at which nanowires with coincid- 
ing sending ends are activated by the C r (one messenger 
particle C r is sufficient to activate nanowires at all the 
six initial positions). Certain ambient gradient can con- 
strain number of the initial positions at which activation 
of the nanowire is allowed (we elucidate that opportunity 
hereafter). 

Temporal aspect of the effective modeling may be ex- 
plained by using a notion of virtual multistate of a com- 



plex molecule that has been activated at a stage T. At 
a step of the effective model process, it is thought as 
appearing in virtual three-state: the state of becoming 
intact by emitting the messenger particle C, the state of 
being probed by sending end of a nanowire at its destina- 
tion position (then, available information about molecule 
activity refers to the stage T only), the state of remain- 
ing intact or becoming active for the future by accepting 
receiving end of the nanowire. Eventually, the molecule 
appears in the single state of being intact or activated at 
the stage T+l. Accordingly, participation of the messen- 
ger particles and the model nanowire in the local infor- 
mation transfer is related only to the virtual three-state 
of the complex molecule as alluded to previously (this 
refers also to virtual bi-state of the complex molecule be- 
ing intact at the T and thus unable to emit the particle 
C). 

The effective model pattern of the local information 
transfer realized by the nanowires can be alternatively 
represented by a formal scheme in which information 
about nanowire orientation, (x, y) in the discrete space 
and the signal Q are coupled into one transferable infor- 
mation portion denoted as Q( x ,y) (see Fig. ^d). Within 
frame of this formal scheme, activation of a nanowire by 
the C r results in appearance of the information portions, 
Q( x ,y) and their distribution to the six target sites being 
the nearest neighbor sites of the x' identified as a site 
nearest the pair [x, y) in the direction from x to y. Thus, 
this x' may be imagined as a node redirecting the hopping 
information portions, Qi x ,y) toward all its closest neigh- 
bors. This imaginary event will be identified hereafter as 
act of hopping- distribution. In accordance with role of the 
signal Q, which has been explained here previously, de- 
livering the Q(x.y) in the way of the hopping-distribution 
to each respective target site allows covering it at a stage 
T + 1 if its nearest neighbor site, so oriented in respect 
to it as x is oriented in respect to y, has been covered at 
the stage T of the evolution. The Q( x ,y) w iH be identified 
hereafter as information portion about covering such tar- 
get site at T+l (we will also use the phrase a signal about 
covering a target site to identify the Q( x ,y))- The scheme 
of the hopping-distribution is used, as more convenient, 
for presenting algorithms of the conditional information 
transfer. 

Let us recall that we model evolution that is random 
process and the scenario presented is to be supplemented 
by introducing certain stochastic condition. Note, in this 
connection, that one of the outcomes of redirecting the 
messenger particle C is reflecting this particle acquiring 
the property C r back to the molecule sending it. One 
may consider a process allowing only this way of redirec- 
tion. Then, only the nanowires whose sending ends coin- 
cide with the site occupied by the molecule sending the 
C can be activated by the messenger particle reflected as 
the C r . This scenario corresponds to the process that can 
be considered as particular case of the process presented 
here previously. In this reduced scenario, the hopping- 
distribution would be the only model mechanism of the 
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FIG. 1: Effective modeiing a path of the focal information transfer in a component system Ai for one step of the evolution, from 
T to T+ 1. (a) The nanowires represented by pairs of neighbor sites are activated by the messenger particle C r and experience 
hopping while being in state of transmitting a signal Q (for the picture clarity, hopping from only one initial position of the 
nanowire is depicted). Three patterns at the bottom illustrate examples of the information transfer model events contributing 
to the virtual effective model event depicted above, (b) Alternative representation of the local information transfer example: the 
Q and the information about nanowire orientation, (x,y) are considered as coupled into one transferable information portion, 
Q(x, y ) ; note that activation of the nanowire does not depend on direction, from which the CV has come (see Section[TIJ. In all 
the patterns: Black circles identify sites covered at stage T; the symbol, ^ detects sites that can be covered for the T + 1 in 
result of the example event depicted; other sites which are indicated by the circles, O are n °t covered at the stage T and will 
not be covered for the T + 1 in result of the example event depicted (those uncovered sites, which coincide with ends of the 
nanowires shown, are indicated by Q). 



information transfer (we call it, direct transfer process 
(DTP) whereas the process encompassing the complete 
distribution of the C r is identified as indirect transfer 
process (ITP)). Note, however, that transfer of the infor- 
mation portion about orientation of a pair of neighbor 
sites in the discrete space is specific feature of the model 
process in both the scenarios, ITP and DTP. For that rea- 
son, one assumes a stochastic condition of occurrence of 
the information transfer realized by the active nanowires: 
in terminology used in the formal model, this is condi- 
tion of occurrence of the hopping-distribution. Moreover, 
assuming the agreement between directions of an active 
nanowire displacements and direction of transmitting the 
Q, from the nanowire sending end toward the receiving 
one, suggests expressing this condition as relationship be- 
tween a value characterizing local, time independent ef- 



fect of stochastic factors and certain ambient gradient of 
a determined forcing factor. Accordingly, we assume a 
constant value parameterizing one determined ambient 
gradient of a forcing factor imposed onto all the ensem- 
ble systems in course of the evolution. Note that using 
a gradient of the forcing factor to express the stochas- 
tic condition will result in constraining a number of the 
initial positions at which activation of a nanowire is al- 
lowed. Value of the gradient pertaining to expression of 
this condition can be assumed so that evolution results 
in stable state without further growth of the lattice area 
covered (see Section IT V AllV|) . 
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III. ALGORITHMS OF CONDITIONAL LOCAL 
TRANSFER 

A. Spatial domain 

All the independent realizations, Ai(T) evolve in the 
same spatial domain under influence of stochastic factors 
and with the same process parameters characterizing the 
model ambient conditions. The spatial domain is a finite 
discrete space x = { x k: k = 1,2, —,N} constituted by 
N < oo nodes of a regular triangular grid embedded into 
a plane so as to fill in a square with the sides of 1 x 1. 
This grid has largest density of nodes among all possi- 
ble grids with the same value of the smallest distance 
between the nodes jl^. Each of the embedded nodes 
is situated at center of a regular hexagon of the honey- 
comb array covering the square. A node as well as the 
hexagon assigned to it are identified also with the term 
"site" (the terms "node" and "site" are used interchange- 
ably here). A node x and its vicinity V x composed of its 
nearest neighbors J/i, J/2, Z/6 are shown in Fig. [2b (the 
V x = {yj , j = 1 -j- 6} together with the x are denoted 
as neighborhood S x ). The A»(T) are represented by pat- 
terns composed only of sites covered at stage T of the 
evolution. 



the stage T, x\ G {Ai{T) n V y ^} (here symbol n denotes 
taking common part of the sets). 

The stochastic condition of the hopping-distribution 
of the Q(x, y ) from the S x is assumed as requirement of 
obeying the inequality P x (y 3 ) > r x (yj,i). Here, r x (yj,i) 
is an element r Xk (^j,i) of the pseudo-random matrix ||r|| 
for x = Xk, where k — 1 -j- N, £j = yj for all j = 1 -f- 6, 
Cj=7 = Xk and i = 1 M. The matrix 1 1 t 1 1 at x y x is 
computed b y u sing a pseudorandom number gene rator 
(RANMAR |3 or, when noted, RANLUX BlU) and 
independently of the stage T of the REP development 
(this is held also for other variants of the processes re- 
ported here). The P x {yj) is the probability assigned to 
a pair at the position (x, yj). The hopping-distribution 
of the Q( x ,y ) which may result from activation of the 
pairs from positions {x,yj) within S x are independent 
events. Here, the value ^ZjP x {yj) < 1 is equal to a 
probability P s of occurring the signal transfers due to 
hopping-distribution of the Q[ x>y .) from the S x and the 
transfer will not take place due to stochastic factors with 
the probability (1 — P s ). 



C. Indirect transfer process REP-ITP 



B. Direct transfer process REP-DTP 

We employ the scheme of the hopping-distribution (see 
Section [HJ to present algorithm of conditional covering 
the sites. Let a node x shown in Fig. [2b be covered in the 
realization Ai at a stage T, x G A{(T) (symbol G is read 
"is an element of). This site is considered active and 
sends signals C that can be reflected only by the neigh- 
bor sites, as signals C r , back to the x. Accordingly to the 
assumption put forward in Section^ a model nanowire 
represented by the pair of neighbor sites at its initial 
position, (x,yj) with yj £ V x , can be activated by the 
signal C r . Index, j identifies here orientation of the pair 
within the S x and pair of indices, jl, with Z = 1,2, 6, 
is used to indicate all target sites, to which the informa- 
tion portions, Q( XtVj ) are distributed if certain condition 
referred to stochastic factors is satisfied. The numbers 
are assigned to the indices in accordance with assign- 
ment shown in Fig. [2b (for clarity, only one example of 
the model nanowire at the initial position (x, yj=z) is de- 
picted in Fig. [2b and reflecting of the C is not shown). 
The denotation, x\ is used to indicate sites being nearest 
the y'ji and so oriented in respect to the y'j l as the site x 
is oriented in respect to yj. In terms used for describing 
scheme of the hopping pairs (see, Fig. QJ,), the denota- 
tion, (x[,y' 3l ) is used to indicate destination positions of 
the pairs. The rule of indexing presented here implies 
the equivalences, x' l=6 = x and y'j l=6 = yj (see the ex- 
ample for j = 3 in Fig. [2b)- The y'^ can be covered for 
the Aj(T + 1) if there is its neighbor site x\ covered at 



This variant of the REP corresponds to a process un- 
derlying evolution of a material system in which infor- 
mation about changing system features at a site can be 
distributed locally to engage more finite size elements 
into the transfer process. 

One assumes that the signals C informing that a site 
z is covered, z G Ai(T) are sent to all closest neigh- 
bors z\ of the z {z\ G V z ) and each z\ redirects the 
signals. In this way, the C is enhanced to C r that is 
delivered to all the x £ V Zl (see Fig. [2b). This communi- 
cation from the z to x is effective if the stochastic condi- 
tion of the hopping-distribution of the Qr x>v A is obeyed, 
Px{yj) > r x(yj,i) and there is the site x\ € {Ai(T)r\V y ', i } 
so oriented in respect to y'^ as the x is oriented in re- 
spect to yj (see Section [ill Bp . In this case, the sig- 
nal, Q{x, Vj ) is delivered to the y'j X that is being covered, 
y'j X G Ai(T+X). This scheme, identified as indirect trans- 
fer process (REP-ITP), allows exploring larger number of 
node configurations close to every site covered at a stage 
T while searching for pairs realizing transfers resulting in 
covering target sites at the stage T + 1 . 

Note that x and z\ can be or can not be elements of 
Ai (T) and this ITP scheme allows also x coinciding with 
the z in the case of reflecting the signal by z\ . Thus local 
realization of the ITP can include the DTP as particular 
case. Eventually, it may happen that x G Ai(T) when 
x 7^ z or x G Ai(T) because x — z. Both the cases 
correspond to contribution of the REP-DTP to the REP- 
ITP. Variant of the REP-ITP that allows only x £ A, L {T) 
will be identified here as pure ITP (REP-ITP-P). 
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FIG. 2: (a) A closest neighborhood S x of a site x £ \ with assignment of the Pj. Formal schemes of the local transfer, the 
example of one finite-size element activated and act of the hopping-distribution is depicted: (b) direct transfer process REP- 
DTP; (c) indirect transfer process REP-ITP; the box with the word, condition and arrow aiming at the big circle means that 
the hopping-distribution shown for example of the pair (1,3/3) occurs if the condition, P x (y3) > r x (y3,i) is satisfied; meaning 
of the circles at sites is explained in the caption to Fig. Q (d) Manifestation of a pattern diversity specific to the REP-ITP 
realization is revealed here by comparing this pattern to a realization pattern at a stage of the spreading percolation SPP (see 
also Section Hvjl . 



IV. PARAMETERS OF AMBIENT 
CONDITIONS 

A. Forcing parameter P 

Ambient gradient of a forcing factor imposed onto all 
systems of the ensemble is one of the macroscopic con- 
ditions. Its effect is modeled in the same way in all the 



realizations Af. by the same specific assignment of the 
probabilities Pj to all pairs (x,yj) within the neighbor- 
hood S x for all the x £ x ( see Fig. E^,-=-c) , 

P 1 = P 6 = P 7 = , P := P 2 = P 3 = Pi = P 5 < 0.25, 

and, accordingly, by initiating the expansion from the 
same pair of sufficiently long vertical chains (we identify 
them as the initiating structure IS; one considers two 
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FIG. 3: Characteristics of the REP advancement shown for examples of REP-ITP representing two classes of evolution: (a) 
without a bend in Ibas(T), here, OL= 0, M = Mq = 100, P = 0.24 (the sum E + 1 means that the values shown should 
be multiplied by ten); (b) with the bend in I bas (T), here, a = |, M = Mo = 100, P = 0.1512. Two of six and two of 
five realizations developing yet at the very late stages are shown for examples of the classes (a) and (b), respectively. The 
realizations are depicted as dark patterns imposed on the brighter pattern of the base. 



chains for the reliability purpose only, see Fig. |2t>) • The 
assignment of the Pj is assumed the same for both the 
REP variants, ITP and DTP, being investigated here. 
The value P parameterizes the model forcing factor af- 
fecting intensity of the expansion (limits to the P varia- 
tion for the REP simulations being studied here are deter- 
mined in Sect ion fVI A|) . The conditions assumed result in 
development of the areas covered within right half-plane 
from the IS (already at first stages of the REP, all the 
realizations evolve from the identical IS into actual finite 
random set, FRS defined in Section HJ. The resulting 
REP tends to stable bi-state represented by the MMS 
patterns that do not grow further (see Sections IV A II 
IV Cf) . The schemes of local covering used to simulate the 
REP result in various forms of the covered area in the 
REP realizations, At. These can be sets of more or less 
dispersed rugged islands (see Fig. [5Ji, |3J) and this diver- 
sity is specific to modeling the finite-size effects. 

This specific feature of the REP realizations is revealed 
by comparison to a realization P>i, with i < M, of the 
spreading percolation (SPP) initiated from the same IS 



and expanding within right half-plane (sec Fig. 01). The 
corresponding assignment of the site covering probabili- 
ties for the SPP is 

Pl = P 6 = , P 7 = 1, P := P 2 = P 3 = Pi = P 5 < 1. 

(2) 

In course of the SPP evolution, a site £ £ S Xi with x E 
P>i(T), is covered for the B t (T + 1) if P x (£) > r x (£, i) 
with £ G {x, yi, 2/2, y&\- Thus, for the SPP, covering 
a site at a stage T implies its covering at T + 1 that is 
not the case for the REP (see Eqs. and Let us 

note that SPP is used here as simplest process, without 
modeling the F-SE and being comparable to the REP, 
which enables us to emphasize effects attributed to the 
F-SE modeling in the REP. 



B. Diversity parameters N and M 

For reliability of the systematic study of the REP, one 
requires comparing results of simulations obtained for ex- 
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tensive series of the parameter values P, in a spatial 
domain x = { x k, k = 1,2, ...,N} with sufficiently large 
number N of elements (nodes) and for appropriate num- 
ber M of the process realizations Aj , i — 1,2, M . 

The number M of sets Ai constituting the ensemble 
(FRS) parameterizes degree to which diversity of the FRS 
manifests itself in features observed in sequence of the 
mean expected patterns, fF(T) (i.e. in the sequence 
of MMS). The value, Mq limiting from below the num- 
ber of realizations M being moderate assures that all 
the realizations have no site covered in common at any 
stage, T > 2 of the random expansion process simulated. 
Number N > Nq fa 18000 of sites constituting the space 
and M > Mq w 100 allow sufficient spatial-specific and 
ensemble-specific diversity modeled by respective com- 
puting the pseudo-random matrices ||r||jvx7xM • The M 
different matrices having the size N x 7 and constituted 
by pseudo-random values are generated, r Xh (£j,i) with 
k = 1,2, ...,7V, ^ = y 3 for j = 1,2,..., 6, £ 7 = x k and 
i = 1,2, ...,M. Assessment of the upper bound Mi to 
the M being moderate results from this investigation and 
is eventually presented in Section IVI B 31 The spatial 
domain, x should be sufficiently large to include each 
realization at latest required stage T mx of the REP de- 
velopment, 

M 

\JS X C X Vx€base[F{T)] = \jAi(T) VT < T mx . 

x i=l 

(3) 

The set, base[F(T)} will be called base of the FRS at the 
stage T of the REP development; the symbols (J and C 
denote union and inclusion of sets whereas V is read " for 
all" . Here, we have found that changes in number N 
only, N <N <16-N do not affect the simulation results 
remarkably. 

C. Lattice order projecting efficiency d 

The finite discrete space \ is an example of a compact 
space that can be covered by finite number of open sets 
having a structural property in common, which is then 
inherited by whole the space (e.g. jl6J. Here, this is 
the form of a regular hexagon. Extensive discussion of 
modelling the lattice order projecting efficiency has been 
presented by us previously |17j . There, we have revealed 
the ways of achieving isotropy of the symmetrical spreed- 
ing on the hexagonal lattice, however, only one of them 
allows regulating smoothly the degree to which the lattice 
order is projected onto a process of the lattice covering. 
This way is followed here: One can select randomly a 
fraction a of sites in x, for which covering all their six 
nearest neighbors is allowed. From the rest of sites, one 
selects randomly a half of them, (1 — a)/2 to allow cov- 
ering only odd sites within their closest vicinities V x (see 
Fig. UK)- For the remained half, one allows covering only 
even sites within the V x . The symmetry and mutual situ- 
ation of the odd and even triplets within the V x result in 



mutual neutralization of their contributions to the inheri- 
tance effect. This can be observed for clusters constituted 
by sufficient number of sites covered. Eventually, the a, 
< a < 1, parameterizes efficiency of projecting the order 
of regular hexagonal array onto patterns resulting from 
the REP development. 



V. CHARACTERISTICS OF THE REP 
DEVELOPMENT 

A. Characteristics of the process advancement 

1. Reference to finiteness and discreteness 

After few initial stages, advancement in the REP devel- 
opment (either ITP or DTP) is characterized by mono- 
tonic growth in the average number of nodes covered in 
a realization, 

M 
i=l 

with Ii(T) — fi{Ai(T)}, and in the number of nodes 
constituting the base, 

Ibas(T) = fi{bas[F(T)]} (5) 

as well as by stepwise growth in the number of nodes of 
the MMS representation, 

I mms {T)=v{jF(T)}. (6) 

Here, denotes the number of elements of the finite 
set {.} (see previous Sections for other denotations). The 
REP finishes at a stage T s , until which all the M real- 
izations have finished their development (see Fig. OK). 
A stage T a , since which I m ms(T) does not grow further, 
identifies beginning of the REP advancement period. Re- 
alizations Ai(T) , i = 1,2, ...,M can finish their devel- 
opment at various stages T < T s . This variety of the 
realization evolutions results in Ibas(T) growing much 
stronger than I av (T) for T>T a (see Fig. |3| . In virtue of 
very small reductions in realizations in all REP simula- 
tions reported here, approaching the T s can be indicated 
by the Ibas{T) and I av (T) becoming constant. 

One may observe specific distribution, h as {T) in the 
beginning of the period following the T a which corre- 
sponds to forms of not overlapping parts of few growing 
realizations. For example, a bend in the has(T), which 
can be observed around a stage T^end occurring close to 
the T a , corresponds to change in form of increments in 
the base from rather uniform front like advancement to 
few narrow separated tongues preserving direction of the 
expansion (see Fig. 12b). 

We have found that smooth transition of the I av (T), 
from strong growth to very weak increments, indicates 
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finishing of the evolution for major part of the realiza- 
tions. These realizations, Aj(T) finish their development 
rather abruptly within a period identified by the series 
of few stages which can be indicated conventionally by a 
central stage, Tb en d{i) of this period. This is indicated 
by a bend in the distribution h(T). The h{T) grows 
strongly until the bend and becomes constant for the 
following stages. Diversification of the values Tb en d(i) 
within the set corresponding to the M realizations and 
lack of the bend in Ii(T) for a number of them result in 
the smooth distribution of the I av (T). 



2. Specific effect of discreteness - steps in I mm .s {T) 

The growing length of steps in the I mm s(T) evidences 
approaching the T a (see Fig.OJ. We explain appearance 
of the steps in the I m ms (T) with reference to the method 
of determining the MMS (see Appendix J. 

At a stage T, one determines the I av (T) and pairs 
of sets, $i(T) or <t>i(T) constituted by such covered 
sites that a number of k > i or k > i realiza- 
tions have them in common, respectively. For the nu- 
merical criterion h = i/M used to determine the MMS, 
one selects such value i that with this i the inequalities 
fj,{(j)i(IT)} < I av (T) < ^{$,(T)} are obeyed. A number 
of such covered sites that more than i — const realiza- 
tions have them in common can grow with T for certain 
period. Therefore, the values /x{$i(T)} and u{(f>i(T)} 
can also grow with the i — const to satisfy the noted in- 
equalities when I av (T) increases with T. At certain stage 
of the following evolution, increment in contribution to 
I av from less overlapping branches of the realizations re- 
quires selecting smaller value i to obey these inequali- 
ties. This change in the i causes jump of increment in 
the Imms as stronger as larger is the contribution from 
not overlapping branches of the realizations. The period 
with i = const is as longer as slower is growth in I av and 
as slower is growth in contribution to the increment in 
I av from dispersed in space, less overlapping branches of 
the realizations. This explains growth in length of the 
subsequent periods corresponding to subsequent values 
i = const, respectively. In virtue of the effects indicated 
by the steps in I mms (T), we will call the corresponding 
time intervals the accumulation periods. 

Thus appearance of the steps results from discreteness 
of the realization ensemble and their significance can be 
attributed to the remarkable share which not overlapping 
extended branches of the realizations have to the I av in 
the REP with the moderate M. 

At earlier stages of the evolution, influence of the var- 
ious factors onto the I m ms (T) is more complex and indi- 
cates significant contribution from overlapping branches 
of the realizations. 



B. Structural characteristics 

Characteristics of local clustering which have been col- 
lected within whole the space x provide information 
about the pattern accomplished within the x at a stage 
T. For that pattern, one may define a probability V n of 
finding a cluster constituted by exactly n sites covered 
within closest neighborhood S x of a site x <E x- We ex- 
press it as V n — ~R.nl N with Mb being number of sites 
within the smallest rectangle circumscribing the base[F] 
and TZ n being number of these sites x with the S x con- 
taining not less than one cluster constituted by exactly 
n sites covered (2 < n < 7); note that using the form 
of a rectangle is just a convenient choice. It is easy to 
see that occurrence of the n-site cluster within S x does 
not allow occurrence of a separate cluster constituted by 
different number, n\ ^ n and ni > 2, of sites within this 

Let us imagine a free gliding window S congruent with 
S x which appears at a site x if the S x coincides with this 
window (see Fig. We follow an approach known 

from the information theory |l8j and consider the win- 
dow S as an information channel. A pattern that can 
be observed within the S can be considered as the in- 
formation carrier passing through this channel. Distinct 
clusters composed of certain number n of sites are inde- 
pendently transfered by the carrier through this channel 
while the free window S appears once at each site within 
the smallest rectangle circumscribing the base[F]. Ac- 
cordingly, one may consider the value, 

H n = -V n ■ InVn (7) 

as contribution from passing the rt-site cluster to mean 
expected amount of information that can be transferred 
by the carrier through the channel S. The sum, H = 
~ Sl=2^« ^ S use d here as reference scale for the H n . 
Distributions of these relationships against the T, not 
just their values at T fixed, characterize the evolution 
of the MMS or Ai realization patterns. They are called 
here, for brevity, the entropic characteristics. We denote 
these values determined for the MMS pattern as, H mm s, 
(H n ) mms , (H n /H)mms- The values, Hi, (H n )i, (H n /H). l 
are determined for a distinct realization Ai whereas H av , 
(H n ) av , (H n /H) av denote mathematical expectations of 
the entropic characteristics determined for every realiza- 
tion of the FRS. 



C. Features of MMS evolution 

1. Structural bi-state 

Accomplishment of the structural bi-state in course of 
a number of initial stages of the REP development is 
specific feature of modeling the finite-size effects. Ap- 
pearance of the bi-state in the MMS development is well 
illustrated by the graphs of the (H n /H) m ms (see 
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FIG. 4: Characteristics of structural bi-state stabilization in the MMS development for example of the REP-ITP with Ot = |, 
M = M = 100, P = 0.1311. The same value T 2 is identified in all graphs (H n /H) mms (T) with 2 < n < 7 (here, only the 
examples are shown). 



Fig. HJ. Saw-shaped distribution of the {H n / H) mms 
against T depicts these values at subsequent stages 
(T, T + l) as corresponding to states of the bi-state. Sit- 
uation of the points representing the subsequent pairs, 
{(H n /H) mms {T) , (H n /H) mms (T+l)} on parallel lines 
for an accumulation period detects structural bi-state 
stability of the MMS development within this period. 
Then, a pair of neighbor stages computed, (T, T + l) 
corresponds to a single physical stage when the system 
may be in one as well as in the other structural state. Re- 
sult of this type has been reported by us previously as 
the one attributed to modeling the finite-size effects and, 
then, its relevance to available experimental observations 
has been noted. 



The REP simulated accordingly to the REP-ITP and 
REP-DTP reveals features of the bi-state stabilization. 
Let us recall that patterns of covered sites, which are 
being accomplished within whole the spatial domain in 
course of the REP-ITP, may be considered as result- 
ing from processes of covering the sites accordingly to 
both the s chemes, REP-DTP and REP-ITP-P (see Sec- 
tion ^^SJ. However, the REP simulations performed 
accordingly to the REP-ITP-P scheme only, resulted in 
appearance of the bi-state with no features of its stabi- 
lization. This shows that sufficient contribution from the 
REP-DTP to the REP-ITP is the necessary condition of 
the bi-state stabilization for the REP-ITP. This condi- 
tion is satisfied for all the REP-ITP simulations reported 
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here (effect of the contribution which REP-DTP has to 
REP-ITP is eiucidated in Section IVTBH . 



REP. Their distributions against P can reveai possible 
regular changes in these evolutions with varying P. 



2. Bi-state stabilization events, To, Ti, T2 

Appearance of the bi-state stability at the end T e of 
certain accumulation period is identified if this effect en- 
compasses the computed sequence from (T e + 1 — 2 • t) 
to the T e , with t > 2, for the (H n /H) mms (T) for all n, 
2 < n < 7. The bi-state stability accomplished is pre- 
served for all the following accumulation periods. One 
may observe such accumulation period that its end in- 
dicated by the T e is the first computed stage in course 
of the evolution, for which the conditions just mentioned 
are obeyed. Then, the stage (T e +1) is denoted as T 2 and 
corresponds to the first jump from one structurally stable 
bi-state to other such bi-state (see Fig. 0J. This corre- 
sponds to a jump in the I m ms{T) between corresponding 
neighbor steps, however, the I m ms(T) cannot replace the 
(H n / H) mms (T) in unambiguous detecting the T2. The 
accumulation period preceding the T2 begins at the stage 
denoted as T\ whereas the one preceding this begins at 
the stage denoted as To. Having given the T2, both the 
Ti and T can be easily detected as stages corresponding 
to local minima in distribution of a helping characteris- 
tic L h (T) = 1 + (H av - H mms )/H av (see Fig. or, for 
higher values of P, as instants of jumps between steps in 

the Imms \F) ■ 

The distributions, (H n /H)i(T) with 1 < i < M, re- 
veal that To or T\ is within a short sequence of the stages 
since which the realizations can start accomplishing the 
bi-state not stable yet or can accomplish stable struc- 
tural bi-state, respectively (see Fig. EJ); hereafter we will 
use also the phrase becoming bi-stable as abbreviation to 
accomplishing the stable bi-state. Significant part of real- 
izations becomes bi-stable until the T2, however, a num- 
ber of them becomes bi-stable very late after the T2 and 
they begin accomplishing the bi-state close to the T\ . 

Eventually, one may consider the triplet, T2, T±, Tq as 
values characterizing evolution of the MMS pattern in the 



J 



3. Characteristic increments in the MMS development 

Sets being increments in the J F(T) for the accumu- 
lation period or at the jump between these periods are 
to be determined for the same state of the bi-state (see 
Fig.|fJJ). These sets are denoted as SjF and identified by 
the pair of indices I, k or Ik for the respective accumula- 
tion period whereas index Tj is used to indicate sets being 
increments at a jump at the Tj. Accordingly, the pairs 01 
or 12 indicate sets being increments in the J F(T) for the 
accumulation periods, from To to T\ or from T\ to T2, re- 
spectively. The accumulation periods, 77. & and instances 
of the jump events Tj are calculated for each state of the 
bi-state. We calculate length of the accumulation period 
between simulated stages for I £ {0, 1} and k = 1 + 1 as, 
n,k = tk—Ti+l, withi fe = Tfe-1 or tk = T k -2, so that t k 
and T; are both odd or they both are even. One considers 
increments SjF^k corresponding to the 77^: between T\ 
and tk for one state and between T; + 1 and tk, A f° r the 
other, with tk, A = tk — 1 if t k =T k — l or t k , A =t k + l 
if t k = T k — 2. The increments S jF^k are determined 
accordingly to Eq. (JSJ. The nonnull reductions, SfF k j 
determined in accordance with Eq. may accompany 
these increments. Changes in forms of these increments 
and reductions with varying the M, Oi and P are consid- 
ered in Section IVI B 31 Increments in JF corresponding 
to a state of the bi-state at the jump in MMS growth at Tj 
are sets of sites scattered throughout whole the J F and 
they are not accompanied by a reduction in the MMS. 
They are denoted as SJ Fy, and determined accordingly 
to Eq. (Unj. Differences A/T(T,T+ 1) between MMS 
corresponding to the two states of the bi-state are given 
by Eq. I|llfl and they have also the form of a set of cov- 
ered sites dispersed throughout the J F. All the sets just 
noted are given by the following relationships, 



6 / F hk e / F(t k )\ / F(T) , / F(t k , A )\ / F(T, + 1) , 



(8) 



SjF k ,i G ljF(Ti)\jF(t k ) , JF(Ti + l)\ jF(t k , A )j, 

6 S Ftj e {/ F{Tj)X I F{Tj ~ 2) ' / F{T] + m / F{Tj ~ i) ] ' 

A J F(T, T+1) = U F(T)\ J F(T + 1) , J F(T + 1)\J T(T)| , 



(9) 



(10) 



(11) 



r 



(here, A \ B denotes elements of a set A which are not MMS depicted in Fig. [H] for example of the 5 JFt 2 as well 
elements of the set P>). The form of increment in the 
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FIG. 5: Examples of characteristics representing three recognized classes of originating the structural bi-state and its stabiliza- 
tion in a realization development. For this example of the REP-ITP with a = |, M = M o = 100, P = 0.1311 (then, T = 14, 
Ti = 20, T2 = 29), of about 35% of the realizations Ai(T) have appeared in stable structural bi-state before the stage T = T2. 
Observe that appearance of the bi-state close to To = 14, at T — 15 > To (i = 75) or at T = 13 < To (i — 43) corresponds 
to its stabilization and becoming the U(T) = const at stages T = 42 > T2 (i = 75) or T — 18 < Ti (i = 43). Here, later 
origination of the bi-state in development of the realization i~57 corresponds to absence of the stabilization in the A57 (T) 
before T = 70 > T 2 . 



as A J F(T, T+l) will be identified as scattered increment. 

Sequences of patterns representing evolution of dis- 
tinct realizations Ai(T) reveal nature of accomplishing 
the bi-state and thus elucidate emerging the characteris- 
tic forms of the increments in MMS in course of the REP 
development. Example of such sequence is presented in 
the array of patterns shown for A 5 j in Fig. Specifi- 
cally developed form of the realization within its central 



region, distanced from the fronts, results from accom- 
plishment of the bi-state. No bias from this specific fea- 
ture in evolution of a realization has been observed cither 
for REP-DTP or REP-ITP simulated here. 
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FIG. 6: Pairs of MMS patterns imposed one on the other to show characteristic increments in the MMS in course of the REP 
evolution. For example of a region in one REP realization the pairs of imposed patterns are shown to reveal appearance 
of the bi-state. In all figures, bright pattern is imposed on the dark one. All figures have been obtained by simulating the 
REP-ITP with a = |, M = M = 100. The top pair of patterns results from the simulation with high P = 0.18 (then, 
T = 82 , Ti=92 , T 2 = 101) whereas all other figures with low P= 0.1311 (then, T = 14 , Ti = 20 , T 2 = 29). 



VI. RESULTS AND DISCUSSION 

One searches for regularities in changes of the MMS 
evolution patterns JF(T) with varying the process pa- 
rameters P, a, M. Results of the previous sec- 
tion suggest that distributions, T (P), T^P), T 2 (P), 
T()i(P), T i2(P) with P varying within an interval when 
a, M are kept constant may reveal such regularities. The 



corresponding increments 5 J Ft^ , 5 J Fqi, 6 J F 12 associate 
the changes in the characteristics with features of the 
MMS evolution patterns. Accordingly, one investigates 
these characteristic distributions with fixed parameters 
a e {0, |, 1} and M = M = 100. Dependence of the 
J F(T) on the number of realizations M is investigated 
with the purpose to estimate the upper limit Mi to the 
values M being moderate. 
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FIG. 7: Regular variations of the MMS bi-state stabilization event characteristics with P and OL. All the REP-ITP simulations 
have been done with M = Mo, the same pseudo random matrix, ||r||Ar x7xm and space of TV = TVo sites. The same pseudo 
random number generator RANMAR has been employed to compute this matrix and to satisfy the condition, OL < 1. Only the 
control distributions shown in the column most to right have been obtained by using the more advanced generator RANLUX 
(with luxury level=3). 



A. Distributions Tj(P), t 01 (P), t 12 {P) 



The three sets of distributions Tj(P), n 2 (P), t 01 (P) 
have been obtained from the series of REP-ITP sim- 
ulations with a £ {1,|,0}, M = M , in the same 
space x {No ~ 18000 nodes) and by employing the same 
matrix 1 1 r 1 1 jv x 7 x Af of pseudo-random numbers. The 
generator RANMAR [p| has been employed to com- 
pute the matrix ||r||jv x7xAf an( i to satisfy condition of 
the lattice order projecting efficiency Ot < 1 (see Sec- 
tion ^^UJ. Three triplets of graphs depicting the distri- 
butions Tj(P) 7 t 12 (P), t 01 (P) with a £ {1, |,0} are pre- 
sented in Fig.0 With the purpose to recognize an effect 
of using a particular pseudo-random number generator, 
they are compared to results of the REP-ITP simulation 
series repeated with employing the more advanced gener- 
ator RANLUX(luxury level=3) [lj, [l5j instead of using 
RANMAR to compute the ||r|| jv x7xM an d to satisfy the 
requirement a = 0. This comparison reveals only quan- 
titative differences that appear to have no effect onto 
inferences about qualitative features of certain regulari- 
ties observed in the distributions. Therefore, the simpler 



generator RANMAR has been considered as sufficient to 
obtain reliable results for the purpose of this research. 

The interval of the values P used to perform series of 
the REP-ITP simulations has been selected specifically 
for each value of the a £ {1, -|,0}. The lower limit of 
the P variations has been selected so as to obtain the 
REP-ITP evolutions revealing clearly occurrence of all 
the three instants T , Ti, T 2 . The upper limit has been 
assumed with the help of the ri2(P). One may observe 
qualitative change in variation of the t\i with growing 
P starting from certain threshold value P = Pt(oi). For 
P > Pt, the value ti 2 (P) appears to vary more regu- 
lar and within fixed interval (see Fig. 0) . This enables 
us to estimate, P T (a= 0) « 0.222 (for RANMAR and 
RANLUX), P T {a= |) w 0.150, P T (a= 1) « 0.060, that 
correspond to the beginning of certain interval [P a ,Pb] 
within which the Tj(P) is growing (see Fig. and Fig. El 
appearance of such intervals is elucidated in next Sec- 
tion ^^^3 ■ The Pt(C() appears to be indicated also by 
specific distributions Tj (P) within two neighbor subinter- 
vals, [P a , Pb) situated just to left of the Pt and within the 
subinterval just to right of the Pr- Eventually, the upper 
limit of the P has been assumed so as to encompass not 
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less than two subintervals [P a , Pb] following the Px{ot). A 
number of test simulations performed accordingly to the 
REP-ITP scheme with the P being substantially higher 
than the limits selected seem to corroborate reliability of 
using the Pt(o) to assess scope of the P abscissa for this 
investigation. 

Form of the graph Tj(P) resembling quasi-periodical 
distribution as well as occurrence of the Pt(cx) are spe- 
cific to the REP-ITP and have not been observed for 
REP-DTP simulation series. The distribution, Tj(P) is 
much less regular for the REP-DTP simulated with the 
M w Mq (see Fig. |HJ and much more complex for the 
REP-DTP series computed with high values of the M. 



B. Structural pattern of changes in REP series 

1. Appearance of intervals [P a ,Pb] 

Interpretation of the distributions Tj (P) requires com- 
paring the corresponding realizations, Ai(T) with the 
same i < M, in the REP simulation series of the same 
type, performed with close values P constituting the 
growing series and with other parameters being the same 
for these simulations. 

One can observe that there are intervals of the value 
P within which a small increment in P may result in 
covering some sites being incipient elements of substan- 
tial branches developing in some realizations. In virtue 
of the M being moderate, those branches can postpone 
occurrence of the T 2 (see also Sections fV A2llV"C 2|) . Fur- 
ther growth in P may result in covering more sites con- 
stituting the incipient elements of the branches in larger 
number of realizations. Evolutions of the corresponding 
realizations obtained with close values of P and for the re- 
spective accumulation periods reveal growth of the spots 
covered, increment in number of those spots and, even- 
tually, spatial uniformization of their arrangement (see 
Fig-IEJ- Then, one may observe growth in number of cov- 
ered spots distributed along an interval being parallel to 
the initial chain (see the example of An in Fig.[HJl or ex- 
pansion of the localized tongues also in this direction (see 
the example of A73 in Fig.lHJ). The spatial uniformization 
of the realizations appears to be sufficient for stabiliza- 
tion of the bi-state at an earlier stage of the evolution 
if P becomes slightly greater than certain value P* (the 
P* appears as right-hand end of the interval alluded to 
previously; the P just mentioned as being slightly greater 
than the P* is denoted as Pi)- This earlier stage appears 
closer to the T\ detected in the evolution performed with 
P being slightly smaller than P* (this value of P is de- 
noted by Pg). Thus the inequality, T 2 (P^) > T 2 (P^) is 
obeyed. Nature of the changes just revealed explains rep- 
etition of the scenario for the REP evolutions with the P 
growing further and thus appearance of subsequent inter- 
vals [P*, P b *} (they are denoted also simpler as [P a , Pb})- 



2. Synchronisation effect of covering sites in REP-ITP 

The distributions Tj{P) obtained from the series of 
the REP-ITP simulations reveal variation resembling 
the quasi-periodical one and T 2 {P) grows rather mono- 
tonically within the intervals [P a ,P,] (see Fig. [7| and 
Fig. . This is accompanied by systematic dependence 
of the average length Cp of intervals [P a ,P,] on the Ct: 
£p(a)/Cp(a— 1) w 3— 2-a. Evolutions of realizations in 
the REP-ITP series simulated with the series of P grow- 
ing within the [P Q ,p,] toward the P, (P — > Pb) result in 
covering larger number of sites dispersed in the space. 
Then, more stages is required to accomplish sufficiently 
dense sets of the sites covered in the realizations to en- 
able sufficient contribution from the REP-DTP to REP- 
ITP for stabilization of the bi-state (this contribution 
has been explained in Section IV C 1(1 . Further growth of 
the P, until Pf,® is reached, results in denser covering 
the area in realizations Aj(T) so that contribution from 
the REP-DTP to REP-ITP becomes more substantial 
and the stabilization occurs earlier. This may be inter- 
preted as result of synchronization between site covering 
accordingly to the REP-ITP-P and REP-DTP schemes 
considered as components of the REP-ITP scheme. 

One may note following factors contributing to the reg- 
ularity of the distributions T.j{P). The REP-ITP differs 
from the REP-DTP in exploring the larger number of 
site configurations within neighborhoods of sites covered 
at a stage T while searching for pairs realizing trans- 
fers resulting in covering target sites at the stage T + 1 
(see Section Ull CI and Fig- Efc)- This results in compa- 
rable contribution of all the realizations {Ai\i—xj-M to 
the MMS in the REP-ITP simulations. Preserving the 
symmetry in the arrangement of the odd or even sites 
whose covering is allowed within vicinities V x with a< 1 
(see Section HV CI and Ref. [IT]]) facilitates more uniform 
covering the space in realizations {Ai}i=i^M and, even- 
tually, appears to improve harmony between the distri- 
butions T 2 (P) , Ti(P) within the intervals [P a ,Pb] (see 
Fig. El - 



3. Variations of MMS evolution patterns with P, OL, M 

Dependence of the increments, SfFoi, SjF\ 2 on the 
simulation parameters P, a with M = Mo is presented 
for the REP-ITP in Table [Q Note that there are no in- 
crements for the accumulation periods following the T\ 2 
if SjFi 2 vanishes. Occurrence of reductions, SjFio and 
S J P21 accompanying increments, S JFqi and 5 jFi 2 in the 
MMS for the accumulation periods toi, t\ 2 as well as 
characteristic differences between those increments (re- 
ductions) in the REP-ITP evolutions, with P — > P a © 
and P — > PbQ for all [P a ,p,], are specific features of 
the REP-ITP simulated with higher values of the Oi; the 
reductions have not been observed for REP-ITP with 
small values of the a and for the REP-DTP (see Ta- 
bles [3[nJ. Appearance of harmony between the distribu- 
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FIG. 8: Spatial uniformization of contributions to the MMS from realizations of the REP simulated with series of the P. This 
uniformization, which corresponds to appearance of the critical value P* , is shown for example of the REP-DTP series with 
M — 104 w M , a= 1 and growing P (here, Pq = 0.10000 < P* < = 0.10125). The figures depict bright pattern imposed 
on the black one. Regular character of distributions Tj(P) within interval of P between subsequent critical values, P a , Pb is 



specific to the REP-ITP; this is shown for example obtained with Oi 



M = Mo = 100. 



tions T (M) , T X (M) , T 2 (M) indicates certain thresh- 
old value M T (see Fig. EH. The REP-ITP simulation se- 
ries done with M slightly smaller or slightly greater than 
M T w 160 result in the specifically different increments 
5 [Fox. Similar effect, concerning however the 8 fFi2, has 
been observed for example of the REP-DTP series (see 
Table Hljl . Let us recall that in course of REP evolution 
simulated with M fss Mo, the increments SjFr appear 
in form of a number of few site clusters and single sites 
scattered throughout whole the J F. On the other hand, 
the corresponding spreading percolation SPP (see Sec- 
tion IIV ^ . being the simplest process without modeling 
the finite-size effects which can be compared to the REP, 
results only in dispersed sedimentation-type increments 
in the MMS that occur only at jumps between steps in 



the I m ms{T) (see Fig. I10[l . These increments can be 
compared to the S J Ft that correspond to jumps in the 
Imms (T) appearing in course of the REP evolution. The 
difference revealed by this comparison can be attributed 
to accomplishment of the bi-state that is peculiar to mod- 
eling the finite-size effects in accordance with the REP 
covering scheme. We have, however, observed that the 
REP-ITP simulated with the sufficiently high M results 
in the S J F^ that, with growing M, are being reduced to 
the forms of dispersed sedimentation-type increments re- 
sembling, to some extent, the increments resulting from 
the SPP (see Fig. I10|l . This resemblance does not mean 
reducing the REP to SPP or vanishing the F-SE; the bi- 
state and its stabilization events at Tj are further present. 
Only effects attributed to discreteness of the informa- 
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FIG. 9: Different degrees of harmony between distributions T2(P) and T\(P) which correspond to two ways of satisfying the 
condition of zero efficiency of the lattice order projecting onto REP (the example of the REP-ITP simulated with M = Mo). 
Figures in the left column depict results obtained with this condition satisfied by random assigning all the possible arrangements 
of sites whose covering is allowed within vicinities V x to all x £ \ . The right column depicts the results obtained with 01= 
satisfied with preserving symmetry in arrangements of sites whose covering is allowed in V x (this method is used here; see 
Section flV (Jl and Ref. \Ll\ for further explanations concerning the CM). 



tion transfer have reduced manifestations in the MMS 
patterns resulting from the REP simulations with the 
sufficiently large M. This reduction affects diversity of 
forms of increments in the MMS. The form of the 6 /-Ft* 
is the one of manifestations of the F-SE which is rele- 
vant to opportunity for employing the REP simulations 
reported here to facilitate designing methods of regulat- 
ing the collective effect accomplishment in corresponding 
real systems (see Section IVI C|) . For this reason only, 
vanishing of this manifestation with growing M for the 
REP-ITP simulation series has been accepted here as cri- 



terion being used to assess upper limit, Mi to the values 
M being moderate. Accordingly, 260 < Mi < 280 has 
been identified by appearance of this qualitative change 
in the increments, SjF^ and SJFt 2 . The distributions 
T 3 (M) detect here also certain value M* w 220 > M T 
such that with, M* < M < Mi, the S J F Tj can be either 
scattered increments specific to the moderate M or may 
resemble the dispersed sedimentation-type increments. 

Let us note, eventually, that results of REP-ITP sim- 
ulations with M « Mi reveal manifestations attributed 
to the discreteness but they are remarkably weaker than 
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FIG. 10: Identifying the upper bound Mi of the M being moderate for example of the REP-ITP series with P = 0.1311, OL = 3. 
For this example, a pair of distributions, l mms (T) illustrates diminishing the effect of discreteness with growing M . Compare 
the increment 8 J Ft 2 (M > Mi) to the increment in MMS at jump in J mm3 appearing in course of the SPP evolution (example 
of the SPP with P = 0.475, OL = |, M = Mo = 100 is depicted by the indicated figure only). The figures depict bright pattern 
imposed on the black one. 



those obtained with much lower moderate M > Mq (see 
Fig. lTU|) . Although, the value of Mi depends on a and P, 
the experience acquired shows that the Mi appears to be 
sufficiently high for performing reliable simulation series 
of the REP-ITP with varying P, a and M < M 1 . Occur- 



rence of the Mi has been detected for the REP-ITP and 
we have not observed the specific changes in 8 J Ft } for the 
REP-DTP simulated even with very high value M=600. 
Inspection of the REP-DTP and REP-ITP realizations in 
course of their evolutions explains this observation: Each 
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TABLE I: Dependence of increments in MMS on P and OL for REP-ITP. In the Tables, Si denotes sedimentation-type increment 
(see Fig. which can be solid or dispersed; symbol identifies a null set or the increment being a single site covered whereas 
~0 shows that the increment is constituted by not more than 2-^4 not clustered covered sites (see Sections IVI IVI Al IVI Bl for 
other denotations). The results presented have been obtained from simulations done with M = Mq. 
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TABLE II: Dependence of increments in MMS on M for exam- 
ple of the REP-ITP with P = 0.1311, a= §, (then M T = 160) 
and REP-DTP with P = 0.104, a= 1 (then M T ss 270 ; see 
Table H] for denotations). 
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the REP-DTP are split into two classes: Each realization 
in the class encompassing of about 70% -7- 80% of all the 
realizations forms only one or two big rugged clusters, 
situated against short subintervals of each IS chain, and 
few very small islands. Realizations constituting the re- 
mained part are represented by very small areas covered 
in their domains and resemble less developed realizations 
of the REP-ITP. Because of this diversification, reduc- 
tion of the increment in the MMS at jumps to the form 
that would resemble this observed for the REP-ITP sim- 
ulated with M ~ Mi , would require values of M much 
higher than the M x detected for the REP-ITP. 
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realization of the REP-ITP is constituted by chains of is- 
lands distributed rather uniformly in direction parallel 
to the initial structure (IS), however, some of realiza- 
tions include narrow elongated tongues (see Fig. . For 
the very high M ps Mi, this results in mean expected 
pattern revealing no changes within central part of the 
J F at the jump. Then, only the frontal parts of the re- 
alizations are sufficiently diversified to contribute to the 
increment in the MMS. On the other hand, realizations of 



C. Discussion and conclusions 

Series of the MMS patterns J F(T) represents accom- 
plishing the model collective effect in course of the REP. 
Discreteness of the REP results in a sequence of jumps 
in the MMS development during the process evolution. 
Modeling of the finite-size effects results in appearance of 
the structural bi-state. Using structural characteristics 
of the MMS pattern enables us to identify stabilization 
of the bi-state at the stage denoted as T2 which corre- 
sponds to certain jump (see Sections IV Bl IV"CT) . This 
T2 is the reference stage used in the search for specific 
changes in features of the MMS evolution pattern with 
varying the process parameters {P, Ct, M}. The tqi and 
T12 separate the triplet of subsequent neighbor jumps in 
values of the MMS characteristics which occur in course 
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of the MMS development at stages Tj G {T , T 1} T 2 ] (see 
Fig. UJ. Occurrence of nonnull sedimentation- type incre- 
ments and reductions in the MMS for the accumulation 
periods toi and t\2 is specific feature attributed to mod- 
eling the finite-size effects in the REP (see Fig. and 
Tables [IJ [HJ . Specific form of the scattered increments 
8 J Ft in the MMS development, which correspond to the 
jumps, is characteristic manifestation of the finite-size ef- 
fects modeling in the REP. Reduction of this particular 
manifestation of the finite-size effects in a series of the 
REP-ITP simulations due to substantial reduction of the 
discreteness manifestations in the MMS patterns, which 
is observed with the sufficiently large M, corresponds to 
the upper bound Mi of the M being moderate. The lower 
bound Mq assures that all the REP realizations have no 
covered site in common and features of the corresponding 
MMS patterns are reliable. 

Accordingly, the MMS evolution pattern JF(T), com- 
puted with the M being moderate Mo < M < Mi, ap- 
pears as sequence of bursts S jFrj in series of patterns 
of the MMS increments. The patterns representing these 
bursts separated by the accumulation periods are con- 
stituted of sites scattered throughout whole the J F(Tj). 
This specific form of the increment <5 jFrj enables us dis- 
tinguishing the d J Ftj from neighbor increments for the 
accumulation periods by performing only observation of 
the MMS increment patterns (see Fig.^J. 

In certain conditions parameterized by, {P, ft, M}, 
a burst S J Ft, separating an accumulation period with 
the sedimentation-type increment in the JF from the 
following one with no or vanishing increment in the J F 
can detect the T\. Here, the examples of such condi- 
tions have been found by performing a number of series 
of the REP simulations (see Tables [U [HJ . These results 
suggest that computing a number of series of the REP 
representing adequately a process underlying evolution of 
certain real system ensemble would enable us to predict 
parameters of real conditions corresponding to values of 
the relevant model parameters, {P, a, M}. In those real 
conditions, stabilization of the structural bi-state of the 
/ F(T) recorded in course of the real process would be 
identified by observing the sequence of the bursts S J Ft, 
and increments in the MMS which may occur for periods 
between them. Note in this connection that identification 
of the increment in the MMS at a jump Tj or for the ac- 
cumulation period refers to the same state of the bi-state 
that can be easily found in the simulation results. In a 
real system, performing this task may require more effort, 
particularly in situation when differences between states 
of the bi-state, A J F(T, T + 1) are comparable to the in- 
crement in MMS at the jump, 8 J Ft . However, for the 
jumps between the stabilized bi-states, the covered area 
representing the difference, A J F(T,T + 1) is distanced 
from the front region of the expanding MMS whereas the 
increment, S J Ft consists of elements occurring within 
whole the MMS (see Fig. Eland Section ES^ . This dif- 
ference between the sets noted is not observed for the 
pure ITP process, REP-ITP-P, in which stabilization of 



the bi-state does not take place. Eventually, the fea- 
ture just recalled would help detecting the T2 in real sys- 
tems corresponding to the REP-ITP or REP-DTP. How- 
ever, the distributions of the Tj(P) with other parame- 
ters fixed are much more regular for the REP-ITP than 
for the REP-DTP with lower M (see Fig. |BJ and simpler 
with higher M. This makes real processes with scheme 
of information transfer corresponding to the REP-ITP 
the first candidates for the regulation. Thus results of 
the computations would enable us reducing a number of 
experiments required to realize the real process with reg- 
ulating the collective effect accomplishment. 

Let us emphasize that the opportunity just revealed 
appears for processes with accomplishing the bi-state en- 
abling us detecting its stabilization events Tj whose role 
has been elucidated here previously. We have also shown 
that realization of transferring the information portions 
by the hopping finite-size elements is a prerequisite for 
appearance of the bi-state. Thus the opportunity afore- 
mentioned will not appear for processes without this kind 
of the information transfer. 

Eventually, we conclude: Pattern of the MMS develop- 
ment, which represents accomplishing the model collec- 
tive effect in course of the REP, experiences predictable 
changes with varying the parameters (P, a, M) specify- 
ing the model ambient conditions. 

This feature is attributed to modeling the finite-size 
effects in the REP simulation series. 

For the REP-ITP, the effective prediction is possible 
with moderate number of realizations, Mq<M <M\. 



D. Final remark 

Let us recall that manifestations of the finite size ef- 
fets (F-SE) appear to be a prerequisite for using changes 
in ambient conditions to regulate accomplishing the col- 
lective effect resulting from evolution of the ensemble of 
discrete systems. The F-SE result from fact that the lo- 
cal information transfer is realized with participation of 
elements represented by pairs of zero dimensional enti- 
ties identified at the nearest sites of the array and expe- 
riencing conditional hopping to positions close by while 
transferring certain information portions and preserving 
their orientation in the space. Diversity of opportunities 
for realizing information transfers by the hopping pairs, 
which is much higher than for only distinct hopping zero 
dimensional particles, required effective modeling whose 
spatiotemporal aspect has been explained in Section iHl 
This has enabled us to reveal the way of regulating the 
collective effect accomplishment. 

On the other hand, it is known that certain evolutions 
of model discrete systems, whose all elements participat- 
ing in the information transfer process allow considering 
them as zero dimensional particles only, may result in 
self-organizing structures composed of sites indicated in 
certain way in course of the process (see eg. Ref. 1}). 
However, the appearing structures do not start acting 
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as the hopping functional elements if no additional as- 
sumptions concerning that functioning in further course 
of the model process are made. This seems to suggest 
that the pairs representing the functional elements are 
required from the beginning of evolution of the ensemble 
systems to afford the opportunity for employing the am- 
bient conditions to regulate the emerging collective effect 
in the way reported here. This particularly concerns the 
task of investigating the influence which changes in ambi- 
ent determined factors can have on collective outcome of 
complex processes taking place in systems of the ensem- 
ble. The possible evolution of swarms of nanoelemcnts 
mentioned here (see Sections |nj as well as turbulent 
transport reported previously by us may be considered 
as the examples of real processes for which this opportu- 
nity for investigating them appears to be useful. 

APPENDIX: MEAN MEASURE SET 
REPRESENTATION 

Below, we present only brief technical way, in which 
the mean measure set representation, MMS is computed; 
the relevant theory is given in |Tlj . 

The MMS , denoted here also as JF has been proposed 
for characterizing mean expected form of a finite ran- 
dom set F — {Ai, i — 1,2, M} constituted by finite 
number M of its realizations Ai , 

jF = {<S> M >h} 1 n{<S>M <h}°, (A.l) 

with the value h £ [0, 1] given by the inequalities, 

1 M 

/i{$m > h} < — 22 »{Ai} < n{$ M > h}, (A.2) 

i=l 

where £ is read "is an element of , denotes the 



number of elements of the set {.}, {<3?m > h} is a set of 
such elements x £ x that ^m(x) > h and symbols {.}°, 
{.j 1 denote set-theoretical cells determined in respect to 
a set {.}. A cell of empty conjunction, D° = {A C 
M : .4 n D = 0} and the cell of inclusion, D 1 = {A C 
Ai : A c n D — 0} are defined as classes of finite sets 
A in the space Ai = 2 X of all subsets of the space %, in 
respect to elements D of the set V covering the space x 
where, V = {Dj C x, 3 = 1, 2, N v } and D k DDj = 
Vj 7^ k. Here, A c denotes complement of the A to M, 
symbol denotes a null set, V is read "for all", R and C 
denote common part and inclusion of sets, respectively. 
The ^>m — &m{x) is the mathematical expectation, 

1 M 

^{x) = —^{x,A i ) l (A.3) 

i=i 

with <J>(x, Ai) = 0if Ai £ {x}° and $(x, Ai) = 1 if A l £ 
{x} 1 , where {a;} is the single element set. Uniqueness 
of the canonical cell representation Eqs. I|A.1|) has been 
proved 01 • The multitude of sets JF is represented 
unambiguously by its single component: if h > 0.5 one 
takes the one whose measure is largest; if h < 0.5 one 
takes the one whose measure is smallest 0]. In this 
work, only such single set representation is used so it is 
denoted in the text by /For MMS. 
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